Nonlinear energy-loss straggling of protons and antiprotons in an electron gas 
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The electronic energy-loss straggling of protons and an- 
tiprotons moving at arbitrary nonrelativistic velocities in a 
homogeneous electron gas are evaluated within a quadratic 
response theory and the random-phase approximation (RPA). 
These results show that at low and intermediate veloci- 
ties quadratic corrections reduce significantly the energy-loss 
straggling of antiprotons, these corrections being, at low- 
velocities, more important than in the evaluation of the stop- 
ping power. 



The energy loss of charged particles has received a 
great deal of attention for several decades fl]-|jj, since 
it plays an important role in investigating the elemen- 
tal composition, depth distribution, and lattice location 
of implanted atoms in matter. The characterization of 
the distribution of energy losses suffered by ions in their 
interaction with matter requires, in the simplest case, 
two quantities: the stopping power and the energy-loss 
straggling. Calculations of the energy-loss straggling of 
charged particles in an electron gas were performed 
within linear response theory, by treating the screened 
interaction to lowest order in the projectile charge (Z\e). 
This procedure leads results for the stopping power and 
the energy- loss straggling that are proportional to Z\e 2 . 
However, early measurements by Barkas et al. on posi- 
tive and negative pions, and also subsequent experiments 
with protons and antiprotons |||, have indicated that 
these quantities exhibit a dependence on the sign of the 
projectile charge. This dependence was ascribed to Z\ 
corrections to the first Born approximation that under- 
lies the linear response theory, assuming that higher odd 
power corrections were negligible in the velocity regime 
under study. Very recent measurements || of the stop- 
ping power for antiprotons in light and heavy targets 
have shown that the stopping power for antiprotons is re- 
duced, near the stopping maximum, significantly (about 
35%) as compared to the corresponding stopping power 
for protons. 

The Z\ correction to the stopping power was first cal- 
culated by using classical perturbation theory for a har- 
monic oscillator |HJ , and quantal calculations of the Z\ 
correction to the stopping power of an electron gas have 



also been performed, in the low |ll]] and high [ p^l3] | ve- 
locity limits. More recently, a quadratic response theory 
for the interaction of charged particles with an electron 
gas has been developed |1J] , in which full account of the 
Fermi motion of the target is taken into account, and 
calculations, within this theory, of the stopping power 
of an electron gas for ions moving with arbitrary non- 
relativistic velocities |l5|] have provided good agreement 
with measurements of the energy loss of protons and an- 
tiprotons in silicon B . A calculation of the Z\ correction 
to the energy-loss straggling and the energy width of the 
states of ions moving in an electron gas has also been re- 
ported, in the low-velocity limit, and it has been demon- 
strated (ljjl^l that this correction to the quantities that 
characterize the distribution of electronic energy losses, 
i.e., the stopping power, the energy-loss straggling and 
the energy width coincides, in this limit, with the result 
of full nonlinear density- functional calculations [pj|-|l9|, 
for high electron densities and small projectile charges. 

In this paper we investigate, within the full random- 
phase approximation (RPA), the Z\ correction to the 
energy-loss straggling of protons and antiprotons mov- 
ing with arbitrary nonrelativistic velocities in an electron 
gas. The electron gas model is well-known to account 
for the contribution to both the stopping power and the 
energy-loss straggling coming from the interaction of the 
projectile with valence electrons in metals po| , and the 
total contribution to these quantities coming from all tar- 
get electrons has also been computed, on the basis of the 
electron gas model of the target, by using the so-called 
local plasma approximation The results throughout 
this paper will be expressed in Hartree atomic units, i.e., 
h = m = e = 1. 

We consider an ion of charge Z\ moving with constant 
velocity v through an isotropic homogeneous electron 
gas embedded in a uniformly distributed positive back- 
ground. The energy-loss straggling per unit path length 
of the projectile, d0 2 /da;, can be expressed in terms of 
the probability of transferring momentum q and energy 
u> to the electron gas, -P(q, ui), as [pOl 



dx v 



w 2 P(q, u))du>. 



(1) 



An explicit expression for this probability is given, up 
to third order in the projectile charge, in Ref. (Ml, by 
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following a perturbation-theoretical analysis of the many- 
body interactions between the moving charge and the 
electron gas. Using this probability one finds, within the 
RPA, from Eq.(f): 

df> 2 fdn 2 \ /dn 2 \ 

where (df2 2 /dx) L and (dil 2 /da;)g represent contribu- 
tions to the energy-loss straggling that are proprotional 
to Z\ and Zf, respectively. (dfi 2 /da;) L represents the 
well-known result one obtains within linear response the- 
ory §§: 

(£) i -4/<^ to <^*»>- < 3 > 

and (d0 2 /dx)g appears as a consequence of the so-called 
quadratic response of the electron gas: 

fdil 2 \ Z\ f d 3 q ~ 0/ . f d 3 q ll; 

x V q _ qi [/i (q,qi)+Mu,Qi) + fa(q,qi)], (4) 

with 

h{q, q x ) = Im^Ree^Ree^ReM , (5) 
f 2 (q, qi ) = Ree-^ee^Ree-i^^g, , (6) 

and 

f s (q, qi ) = -2Im e - 1 Ime- 1 Ree-i giJ ff| i?1 . (7) 

Here, q = (q, u>), q\ — (qi,wi), V q is the Fourier trans- 
formation of the electron-electron bare Coulomb interac- 
tion, 



t q is the longitudinal dielectric function of the medium, 
M* represents the quadratic response function, H* 
is related to the imaginary part of M* , as shown in 
Rcf. Jlij] , 6(x) is the Heaviside function, ui — q • v, and 
wi = qi • v. 

We have calculated linear and quadratic contributions 
to the electronic energy-loss straggling for a wide range 
of non-relativistic velocities of the projectile, after sub- 
stitution of the full RPA linear and quadratic response 
functions into Eqs. (3) and (4). The result of this 
calculation for an electron density parameter r s = 2 
(r s = [3/47T7J] 1 / 3 , n being the electron density) and 
Z\ = 1 is presented in Fig. 1, as a function of the veloc- 
ity of the projectile, by dashed (linear contribution) and 
solid (quadratic contribution) lines, together with the to- 
tal (linear + quadratic) energy-loss straggling of protons 
(dashed-dotted line) and antiprotons (dotted line). It is 



obvious from this figure that at low and intermediate ve- 
locities quadratic contributions result in a significant re- 
duction of the energy- loss straggling of antiprotons. On 
the other hand, though at low velocities the Z\ effect 
appears to be more important in the energy-loss strag- 
gling than in the evaluation of both the stopping power 
and the energy-width, at high velocities the quadratic 
contribution to the energy-loss straggling decreases very 
rapidly. 

The numerical results of the total RPA Z\ contribu- 
tion to the energy-loss straggling for different values of 
the electron density parameter r s are illustrated in Fig. 
2. At low velocities, the Z\ term gets larger as the elec- 
tron density increases, as discussed in Ref. [jl6|. At high 
velocities, Z\ corrections to both the stopping power and 
the energy-loss straggling present similar dependences on 
the electron density. 

For completeness, the low-velocity limit of the energy- 
loss straggling of protons and antiprotons is presented in 
Fig. 3, as a function of the electron density parameter r s , 
as obtained up to third order in the projectile charge af- 
ter substitution of the low-frequency forms of both linear 
and quadratic response functions into Eqs. (3) and (4). 
Since exchange and correlation are absent in our RPA 
treatment, we have also represented, in the same figure, 
full nonlinear self-consistent Hartree calculations like the 
ones presented in Ref. This figure indicates that 

in the case of antiprotons with r s ~ 2 our quadratic re- 
sponse calculations agree, in the low velocity limit, with 
the corresponding full nonlinear results; accordingly, our 
calculations are expected to give, for these electron densi- 
ties, good account of the full nonlinear energy-loss strag- 
gling of antiprotons moving at arbitrary velocities. In 
the case of protons our quadratic response calculations 
appear to be accurate, in the low velocity limit, only at 
high electron densities (r s ~ 1.5). 

The results depicted in Figs. 1, 2 and 3 indicate that 
the Z\ contribution to the energy-loss straggling in an 
electron gas is positive for Z\ > 0, for all velocities of 
the projectile and all electron densities of the target. As 
a result, the energy-loss straggling is, within quadratic 
response theory, greater for a proton than for an antipro- 
ton, as in the case of the stopping power. This result is, at 
high electron densities, in agreement with full nonlinear 
self-consistent Hartree calculations, as discussed in Ref. 
|]l6| . At low electron densities self-consistent Hartree cal- 
culations still predict, in the low- velocity limit, nonlinear 
effects to reduce the energy-loss straggling of antiprotons. 
However, the formation around protons of bound states, 
not described within perturbation theory, tends to screen 
out interactions with the electron gas, and this causes, 
at low electron densities, a significant reduction in the 
energy-loss straggling of protons, as can be inferred from 
Fig. 3. This reduction can be accounted for through the 
introduction of an effective charge, as discussed in Refs. 
@and@. 

In summary, we have calculated Z\ corrections to the 
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energy-loss straggling of protons and antiprotons moving 
at arbitrary nonrclativistic velocities in a homogeneous 
electron gas, by using a quadratic response theory and 
the RPA. At high velocities (after the plasmon thresh- 
old), the quadratic response theory provides an accurate 
estimate of the full nonlinear correction to the energy- 
loss straggling of protons and anitprotons moving in an 
electron gas. In the case of antiprotons moving at arbi- 
trary nonrelativistic velocities in many solids used in the 

experiments (r s ~ 2), the present quadratic response the- 
ory provides an accurate description of nonlinearities in 
the energy-loss straggling. A comparison with results ob- 
tained within linear response theory indicates that at low 
and intermediate velocities the Z\ correction reduces sig- 
nificantly the energy-loss straggling of antiprotons, this 
correction being, at low velocities, more important than 
in the evaluation of the stopping power. 
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FIG 1. Linear (dashed line) and quadratic (solid line) con- 
tributions to the energy-loss straggling per unit path length 
for Zi = 1 and r s = 2, divided by the velocity of the projectile 
and as a function of the velocity. Total (linear + quadratic) 
energy-loss stragglings of protons and antiprotons, divided 
by the velocity, are represented by dashed-dotted and dotted 
lines, respectively. 

FIG. 2. Quadratic (Zf) contribution to the energy-loss 
straggling per unit path length, divided by the velocity of 
the projectile and as a function of the velocity, for Z\ = 1 
and four representative values of r s : r s = 1, r s = 2, r s — 4, 
and r s — 6. 

FIG. 3. Total (linear + quadratic) energy-loss straggling 
per unit path length divided by the square of the projec- 
tile velocity, as a function of the electron density parameter 
r 3 , for protons (dashed line) and antiprotons (dashed-dotted 
line). Crosses and stars represent the full nonlinear en- 
ergy-loss straggling of protons and antiprotons, respectively. 
The solid curve represents the linear contribution to the en- 
ergy-loss straggling for Z\ = 1. 
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